Geometric phases and competing orders in two dimensions 



Liang Fu and Subir Sachdev 
Department of Physics, Harvard University, Cambridge MA 02138 

Cenke Xu 

Department of Physics, University of California, Santa Barbara CA 93106 

(Dated: October 20, 2010) 



O 

(N 
-t— > 
O 

O 

00 



u 

■i— > 

B 

local antiferromagnetic order, induced by geometric phases. We also discuss the connection to 
topological terms obtained under situations where the electronic spectrum has a Dirac form. 
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Abstract 

We discuss the problem of characterizing "quantum disordered" ground states, obtained upon loss 
of antiferromagnetic order on general lattices in two spatial dimensions, with arbitrary electronic 
band structure. A key result is the response in electron bilinears to the skyrmion density in the 
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I. INTRODUCTION 



Geometric phases have played a central role in many fields of physics, and especially 
in the quantum Hall effect at high magnetic fields^. However, in subsequent research it 
has also become clear that geometric phases are crucial for a complete understanding of 
the quantum phase transitions of correlated electron systems in zero applied magnetic field. 
Traditionally, classical phase transitions are described in terms of an 'order parameter', with 
one phase being ordered, and the other 'disordered'. Upon extending this idea to quantum 
phase transitions, we have the possibility of a 'quantum-disordered' phased However, in 
almost all of the interesting examples, the latter phase is not disordered: geometric phases 
induce a 'competing' order. A separate possibility is that the quantum-disordered phase 
has fractionalization and topological order: we will not explore this latter possibility in the 
present paper. 

(We note here that the word 'phase' has two separate meanings above, and in the remain- 
der of the paper. When used by itself, 'phase' refers to a particular state of a thermodynamic 
system. However, in the combination 'geometric phase', it refers to the angular co-ordinate 
of a complex number representing the wavefunction. We trust the context will clarify the 
meaning for the reader.) 

In two dimensional systems, the earliest example of a competing order induced by geo- 
metric phases was in the spin S = 1/2 square lattice Heisenberg antiferromagnet. The model 
with nearest-neighbor interactions has long-range Neel order. We can try to destroy this 
order by adding further neighbor frustrating interactions, leading to a possible quantum- 
disordered phased Such a phase should be characterized by the proliferation of defects in 
the Neel order: for collinear Neel order in a model with SU(2) spin symmetry, the order 
parameter lies on S 2 (the surface of a sphere), and the homotopy group ^(S 2 ) = Z, implies 
that the existence of point defects known as hedgehogs. Haldane^ pointed out that geometric 
phases of the lattice spins endowed each hedgehog with a net geometric phase, and argued 
that this implied a 4-fold degeneracy of the quantum-disordered ground state. Read and 
Sachdev^El demonstrated that the hedgehog geometric phase actually implied a competing 
order, associated with a broken lattice symmetry due to valence bond solid (VBS) order. 
The VBS order can take 4 orientations related by lattice symmetries, thus realizing the 4- 
fold degeneracy. They also presented two additional derivations of the hedgehog geometric 
phase: from the Schwinger boson representation of the spins 5 , and via a duality transform 
of the quantum dimer modeP. Sachdev and JalabertP introduced a lattice gauge theory for 
the competing Neel and VBS orders, in which the geometric phase appeared as a coupling 
between the skyrmion density associated with the Neel order (defined in Section [H]) and a 
lattice field linked to the VBS order. 

We note in passing that we will not be interested here on the separate question of the 
nature of the phase transition between the two competing order phases. A second order 
transition appears in the 'deconfined criticality' theory^, and this proposal has been the 
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focus of a number of numerical studies^^. 

A different perspective on the Neel-VBS transition appeared in the work of Tanaka and 
HrP^ who used a continuum theory of Dirac fermions. The previous works had all represented 
the spins in terms of bosonic degrees of freedom which carried geometric phases. Tanaka 
and Hu instead used a fermionic representation of the spins, and then chose a low-energy 
limit which allowed representation of the theory in terms of continuum Dirac fermions in 
2+1 spacetime dimensions. The Dirac representation appeared from a band structure of the 
lattice fermions in which there was tt flux per plaquette: this could be interpreted as the 
dispersion of fermionic spinons in a a particular algebraic spin liquid (ASL) important for 
intermediate length-scale physics, or as a mean-field dispersion of electrons in a particular 
extended Hubbard modeP^. Armed with the Dirac fermions, Tanaka and Hu used field- 
theoretic developments by Abanov and Wiegmanri^ to show that the effective action for the 
Neel and VBS order parameters allowed representation of the geometric phase as a Wess- 
Zumino-Witten (WZW) term for a 5-component order parameter. The co-efficient of this 
WZW term was quantized to a value reduced consistently to the hedgehog Berry phases 
in the appropriate limit. In a different context, Grover and SenthiP^ recently showed that 
a WZW term was also present between a quantum spin Hall order parameter and s-wave 
superconducting order on the honeycomb lattice; their computation also used the Dirac 
spectrum of the electrons on the honeycomb lattice. 

The appearance of the WZW term with quantized co-efficient in the above computation 
appears surprising from the perspective of the earlier bosonic formulations^'. In these earlier 
works, the quantization was directly related to the quantization of the spin on each lattice 
site, which relied crucially on the projection to one electron (in general, to 2S electrons 
in a fully antisymmetric orbital state, for spin S) on every site. In contrast, in the above 
fermionic formulation ^ 18 1 20 1221 , the local constraints are ignored in the computation of the 
WZW term, apart from a global constraint on the average fermion density. We will argue in 
this paper that the WZW term with a quantized co-efficient is an artifact of the low energy 
Dirac fermionic spectrum. 

This aim of our paper is extend the use of the fermionic representation of the underlying 
degrees of freedom to cases without a low energy Dirac limit. We will develop a general 
approach to computing geometric phases, which works for arbitrary electronic band struc- 
tures, whether insulating, metallic or superconducting. Like the recent workp2GSG2 we will 
not impose a local constraint on the electron number, which is permissible for the metallic or 
superconducting cases or even in insulators with small on-site repulsive energy. Our compu- 
tation begins by applying local antiferromagnetic order, and computing the band structure 
in the presence of this order. Then, we allow the orientation of the local order to become 
spacetime-dependent, so that eventually there is no true long-range antiferromagnetic order. 
However, the local ordering is still assumed to be present, with its associated band structure, 
and we fill these electronic states up to the Fermi level. We will then compute response of 
these filled electronic states to spatial variations in the antiferromagnetic order. We will also 
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allow spatial variations in competing orders, deduce their coupling to antiferromagnetism. 
We will find geometric phases between the order parameters, but will show that a WZW 
representation does not exist in general. 

Another approach to the general problem of geometric phases was described recently 
by Yao and LeeP Their method required extension^ of the 2 dimensional electronic band 
structure to 6 dimensions, and the computation of topological invariants in 6 dimensions and 
of the mapping between 2 and 6 dimensions. Non-zero values of these invariants were then 
argued to be sufficient conditions for a WZW term in the effective action for the competing 
orders. This last conclusion appears to be at variance with our results. 

We begin in Section [TT] by considering spatial variations in the antiferromagnetic order 
on the square lattice. We compute the response to this spatially varying background, in the 
spirit of the computation of Chern numbers of integer quantum Hall states by Thouless et 
a/P. This leads to the key result in Eq. d2.22|). 



Section III extends the computation to allow for simultaneous variation of both Neel and 



VBS orders. Here we will also make a connection to the dimensional reduction methocP^] 



noted above. Section IV contains applications of our results to insulators on the honeycomb 
lattice, while Section [V] considers transitions in the background of the nodal quasiparticles 
of a (i-wave superconductor. 

II. FLUCTUATING NEEL ORDER 

Our approach begins with with an arbitrary band structure for lattice fermions c a , with 
the spin index a =t, U so the band structure of the electronic quasiparticles is 

ff 6 = -J]t(r i -r i ) C t(r < )c a (r i ) (2.1) 

where labels the lattice sites, and t(r) are the tight-binding hopping matrix elements. For 
defmiteness, let us consider the Neel state on the square lattice, as described by the Slater 
mean-field theory of antiferromagnetic order. We allow the Neel order to have a slow spatial 
variation in its orientation, which we specify by the unit vector n a (r) (a = x,y,z). In this 
modulated Neel state, the electronic quasiparticle Hamiltonian is modified from the band 



structure in Eq. (2.1) to 



H = -J2^ rt ~ T i) c U T M T i) +mJ2vin a (r i )cl(r i )a^(r i ) (2.2) 



where a a are the spin Pauli matrices, ^ = ±1 on the two sublattices of the Neel order, and 
m is a mean-field magnitude of the band splitting due to the Neel order. The main result 



of the following Section II A will be obtained by working directly with Eq. (2.2) for a slow 



variation of n a (r) about a fully polarized Neel state. 
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(2.3) 



For some purposes, we will find it advantageous to use an alternative gauge-theoretic 
formulation, which has some technical advantages for a global perspective on the phase 
diagram. For this, we follow Ref. |2BJ and transform to a rotating reference frame in the 
varying Neel background so that the Neel order points in the constant direction (0,0,1) 
in the new reference frame. We do this by introducing complex bosonic spinors Zi a , with 
\z^\ 2 + \zii\ 2 = 1 so that 

/ Ct | = ( Zt 

where ip p , p — ±, are the "electrons" in the rotating reference frame. We will assume that 
the z a have a slow dependence upon spacetime, allowing in expansion in gradients of the z a . 
A fixed orientation of the Neel order is realized in the rotating reference frame by choosing 
the z a so that 

n a = z* a o%zp (2.4) 

However, we will not assume any slow variations in the fermions c a and i(j p , allowing them 
to carry arbitrary momenta and band structures. 




Parameterizations like (2.3) were motivated earlier by the Schwinger boson formulation of 
the underlying antiferromagnet. In such theories, the geometric phases of the spins at half- 
filling were associated entirely with those of the Schwinger bosons^. In our computations of 
geometric phases in the present paper, we will find it convenient to work in an approach in 
which the lattice geometric phases are attached entirely to fermionic degrees of freedom. For 
this, we will use an exact rotor model formulation of a general lattice Hamiltonian for which 



Eq. (2.3) also holds. The details of this rotor formulation are presented in Appendix |AJ and 
this should be regarded as an alternative to earlier Schwinger boson formulations. In the 
rotor theory, the ip± are canonical fermions with a density equal to the full electron density; 
thus in the the insulator, the total if>± density is 1, and it is this unit density which leads 
to the geometric phases. The bosonic variables z a have a rotor kinetic energy with only a 
second-order time-derivative in the action i.e. they are not canonical bosons, and do not 
directly carry any geometric phases. In the Schwinger boson formulation, the bosons are 
canonical, and this complicates the computation of geometric phases in the general case. 



Inserting Eq. (2.3) into Eq. (2.2), we obtain the theory for the ip± fermions, which we 



write in the formpS 

H = -^tfa ~ *M^y pAij ^ r i) +mY t ViP^ P (r i )M T i) + ■■■ ( 2 - 5 ) 



First, note that the transformation to the rotating reference frame has removed the slowly 
varying r dependence from the second term proportional to m. Instead the effect of the 
transformation into the rotating reference is now entirely in the hopping term. As discussed 
in earlier workP^, these modifications can be expressed in general in terms of a SU(2) gauge 
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potential, corresponding to the SU(2) gauge redundancy introduced by the parameterization 



in Eq. (2.3). In the fluctuating Neel state we consider here, the SU(2) gauge invariance 



is 'Higgsed' down to U(l): this corresponds to the invariance of Eq. (2.4) only under a 
U(l) gauge transformation of the z a . So we write only the U(l) gauge potential term in 



Eq. (2.5), represented by Aij. The ellipses in Eq. (2.5) refer to additional fermion hopping 



terms connected to the remaining SU(2) gauge fields: these were written out explicitly in 



Refs. I2"6"f27j and also appear in the present paper as the last two terms in Eq. (A21). 



As we are using a continuum formulation for the order parameter n a (r) and the z a , we 
should also work with a continuum U(l) gauge potential A(r). This is related to by an 
integral on straight line between and 



= duA (i-j + u{vj - Yi)) ■ {vj - r % 
Jo 



(2.6) 



The flux in the continuum gauge field A can be related to the 'skyrmion density' in the 
antiferromagnetic order parameter: 



d x A y 



d y A x 



e abc n a d x n b d y n c . 



(2.7) 



With periodic boundary conditions, the spatial integral of the skyrmion density on the right- 
hand-side is a topological invariant, and is quantized to an integer multiple of 2ir; the integer 
is the skyrmion number. Thus inducing a 2tt flux in A corresponds to changing the skyrmion 
number of the field n a (r) by unity, which is the same as introducing a hedgehog defect in 
the Neel order. 



A. Response to spin textures 

This section will carry out the formally simple exercise of computing the linear response 



of the Hamiltonian Eq. (2.2) a slowly varying spacetime dependence in the order parameter 
n a (r). A similar computation can also be carried out using the alternative gauge-theoretic 
form in Eq. (2.5) to a slowly varying gauge potential A^: the latter computation is presented 



in Appendix [B] 



We begin with the Hamiltonian in Eq. (2.2), and assume n a (r) is a slowly varying unit 



vector. In any local region, without loss of generality, we can choose co-ordinates so that 
n a (r) is close to the North pole (0,0,1). In this co-ordinate system, as in Ref. [27, we 
parameterize the variations in the Neel order in terms of the complex field tp via 



n 



if + (p* if - if* 



2i 



m 



(2.8) 



We assume \<f\ <C 1 and slowly varying. Inserting Eq. (2.8) into Eq. (2.2) we obtain the 
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Hamiltonian H = Ho + Hi with 

Ho = J2 (^c f (k)c(k) + mc f (k + Q> 2 c(k)) , (2.9) 

k 

where Q = (ir,ir) and 

e k = -^t(s)cos(k-s), (2.10) 

s 

with t(— s) = t(s). Throughout this section, the summation over momenta extends over the 
entire square lattice Brillouin zone. Also, we will drop the a spin indices of the c a , all Pauli 
matrices in this present section will be assumed to act on the a space, and the a indices will 
be traced over. The coupling to the spatial variations in the Neel order parameterized by tp 
are given to the needed order in (p by 

H x = m [^(kOc^ka + Q>+c(k 2 + kx) + ^(kOc^ka + Q>- C (k 2 - k t ) 

ki,k2 

-y ^(ki)v?(k 1 + k 2 )c t (k 3 + Q)a 2 c(k3-k 1 ) (2.11) 

ki,k2,k3 

We are now interesting in computing the response of the observable properties of H to a 
slow variation in the Neel order n a (r). A key choice we have to make here is that of a suitable 
observable. We are interested in the nature of the phase where Neel order is 'disordered' 
and so it is natural that the observable should be spin rotation invariant. Also, because we 
will use the observable to characterize a 'competing order', it should preferably vanish in 
the spatially uniform Neel state, and be induced only when there are spatial variations in 
the Neel order. Finally, for convenience, the observable should be a fermion bilinear. With 
these constraints, it turns out that a unique choice is forced upon us: it is the observable 

0(k, r) = I <c f (k + Q + q/2)c(k - q/2)) e^ qr (2.12) 

Here the integral over q is over small momenta, characteristic of those carried by the bosonic 
fields; thus the variation of C?(k, r) with r is slow. In the simplest case, the right- hand- side 
has support only at q = 0, so that (D(k, r) takes the r-independent value 

0(k) = (ct(k + Q)c(k)) . (2.13) 

On the other hand, k is an arbitrary momentum in the Brillouin zone, and we will find very 
useful information in the k dependence of 0(k). It is easy to check from H that C?(k) = 
in the uniform Neel state, as we required; only ^d"(k + Q)cr 2 c(k)) ^ in the uniform Neel 
state. We present an alternative derivation of the choice of the observable O in Appendix |B| 
there we consider an arbitrary fermion bilinear, and show that it is O which is uniquely 
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induced to leading order in the applied gauge flux. 

We now proceed to a computation of 0(k, r) in powers of ip using the Hamiltonian 
H + Hi. We will need to work to second order in (p, and also to second order in spatial 
gradients of (p; as stated earlier, all fermion momenta are allowed to be arbitrary at all 
stages. 

First, let us collect the propagators of Ho. The single fermion Green's function of Hq is 
written in terms of its 'normal' and 'anomalous' parts as 

(c(k); C t(p)) = 5 kiP G(k) + 5 k+Q , p a z F(k) 
G(k) = —^r- + 



iuj + Eik —iw + E 2 \ L 

F(k) = u k v k ( , 1 ■ * ), (2.11) 

V —iu + Ei k -iuj + E 2k 



where k takes all values in the square lattice Brillouin zone. The eigenenergies in Eq. (2.14) 

are 



£i, 2k = 2 " ± y 2 - /» 2 - (2.10) 
and the parameters are 

w k = cos(# k /2) , v k = sin(fl k /2) (2.16) 

with 

TYl 

tanfl k = - — , 0<^ k <7r (2.17) 

Note that these relations imply 

Mk+Q = ^k , ^k+Q = «k , -&L,k+Q = E lk , £ , 2,k+Q = -^2k- (2-18) 

The contributions to (c^(k + Q + q/2)c(k — q/2)) to second order in <p are shown in Fig.jl} 
The last diagram vanishes identically, while the first two evaluate to 

( C t(k + Q)c(k + qi - q 2 )> = J ( k ' <h)<P*{<toM<li) (2-19) 

qi,Q2 

where 



J(k,qi,q 2 ) = 



F(k)G(k + Q-q 2 )G(k + qi -q 2 ) 



G(k + Q)F(k - q 2 )G(k + qi - q 2 ) + G(k + Q)G(k - q 2 )F(k + q x - q 2 ) 



- F(k)F(k-q 2 )F(k + qi -q 2 ) 



- (qi -q 2 ) (2.20) 



S 




FIG. 1: Diagrammatic perturbation theory for O using the couplings in Hi in Eq. (2.11). The 



wavy lines are ip sources, the filled circle is the O source, while the full lines are c propagators. 



We now expand this to second order in q x and q2. This leads to very lengthy expressions, 
which we simplified using Mathematica. In the end, a simple final result was obtained: 



J(k,qi,q 2 ) = (qx X q 2 ) 



dk 



dk 



m 



iu + E lk ) 3 (-iuj + E 2k ) 



(2.21) 



Now we combine Eqs. (2.19) and (2.21). The Fourier transform of (qi x q2)<^*(q2)y?(qi) is 



d x (fd y ip* — d y ipd x ip* and to second order in ip this equals —2i(d x n x d y n y — d x n v d y n x ). In a 
spin rotationally invariant form, this expression is proportional to the skyrmion density, and 
so we have one of our main results: 



C(k,r) = -iF{k)e abc n a {T)d x n\v)d y n c {v) 



(2.22) 



where 



T-(k) 



X 



dk 



dk 



2m 3 



-ico + £i k ) 3 (-2w + E 2 k) 



3 I'dek+Q de k \ (sgn(£ lk ) - sgn(£ 2 k)) 



(2.23) 



dk dkj (E lk -E 2k )5 
In the last step, we have evaluated frequency summation at zero temperature. In the re- 



9 



maining analysis we will assume we are dealing with a fully gapped insulator with > 
and i?2k < over the entire Brillouin zone. The metallic case has singularities at the Fermi 
surfaces which are at Eiy_ = or = 0, but we will not explore its consequences here; in- 
deed in our expansion in powers of qi^, we have implicitly assumed smooth behavior across 
the Brillouin zone. Note that in both the insulator and the metal there is no singularity due 



to the denomination in Eq. (2.23): via Eq. (2.15) we always have Ei^ — i?2k > 2m. 



A plot of .F(k) for the insulating case is shown in Fig. U, 




FIG. 2: A plot of the function -^(k) in Eq. (Bll) for = cosk x — cosk y + 0Acos(k x + k y 
0Acos(k x — k y ) and m = 1. 



+ 



The integral of ^(k) is zero over the Brillouin zone. However, note that it has the 
same symmetry as the function (cos k x — cos k y ) sin k x sin k y ; so the integral of ^(k) (cos k x — 
cos k y ) sin k x sin k y is non-zero. This suggest we define the charge Q by 



<2 = — i c^(k)c(k + Q)(cos k x — cos k y ) sin k x sin k y . 



(2.24) 



Note Qt = Q, 



Our main result in Eq. (2.22) implies that any quantum fluctuation which leads to a non- 



zero value of the skyrmion density e a b c n a d x n h d y n c will induce a change in O. Generically, 
a change in O must imply a corresponding change in Q because the two observables have 
identical signatures under all symmetries of the Hamiltonian. In paticular, a hedgehog 
tunneling event is one in which the spatial integral of e a ^ c n a d x n b d y n c (the skyrmion number) 
changes by 47r. Thus, before the hedgehog event (Q) = 0, while after the hedgehog tunneling 
event, we have (Q) ^ 0. We can normalize Q so that (Q) = 1 for each hedgehog, and the 



normalization constant will depend upon Eq. (2.23) and the details on the band structure. 
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Then with such a normalization, we have the important correspondence 

Q = skyrmion number. (2.25) 

This is the key result of the present subsection. We emphasize that such a correspondence is 
possible because both the skyrmion number and Q are invariant under spin rotations, have 
identical transformations under all square lattice space group operations, and are both odd 
under time-reversal. 



B. Connection to VBS order 



The results in Eqs. (2.22) and (2.25) suggest strong consequences in the 'quantum dis- 



ordered' phase where Neel order has been lost. Such a phase will have a proliferation of 



hedgehog tunnelling events, and so Eq. (2.25) implies that there will be correspondingly 



large fluctuations in the charge Q. We can therefore expect that fluctuations in variables 
conjugate to Q will be suppressed, and will therefore have long-range order: this is the 



competing order induced by the geometric phase in Eq. (BIO). Thus any quantum variable 



conjugate to Q is a bona-fide competing order. There are many possibilities, but here, we 
verify that the traditional VBS order does satisfy the requirements. A more specific field- 
theoretic discussion of the appearance of VBS order in the quantum-disordered Neel phase 



will be given in Section III A 

The VBS order is V — V x + iV v defined by 



V x = i c f (k)c(k + Q 8 ) sin k x 

k 

V y = zy^c t (k)c(k + Q y )sin k y 



(2.26) 



where Q x = (vr, 0) and Q y = (0,7r). Now we can compute the commutators 

-m +/, \ /, ^ \ . , (cos(k x ) — cos(3k x )) 
[Q,V X ] = -JJc^kMk + Q^suiV — — - lV v 

k 

[Q, V y ] = ct ( k M k + Q») sin fc* {C ° S{ky) ~ C ° S(3ky)) * -iV x (2.27) 



Here the ~ means that the two operators have the same symmetry under the square lattice 
space group. Thus we have the key result 



[Q, V] ~ V. 



(2.28) 



This means that V is a raising order for Q. But, as we noted in Section II A this is precisely 
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the effect of the monopole tunneling event: in other words, V has the same quantum numbers 
as a monopole operator. Then, following the reasoning in Refs. 5 9, we may conclude that 
V is a competing order which becomes long-range in the quantum-disordered Neel phase. 
An alternative route to determining an operator conjugate to Q is to determine a V so 



that —i{y^d t V — Vd t V^) ~ Q. It is easy to check that the definition in Eq. (2.26) does 
satisfy the needed requirements. We have the time derivative 



dt 



E 

k 



sin fc K (£ k - £ k+Q Jc f (k)c(k + Q x ) + 2m sin k x c j {k)a z c{k + Q y ) (2.29) 



and similarly for V y . For simplicitly, we will drop the terms proportional to m, and work in 
the limit of small m. So we have 



~ ' '( Vy-jr - Vx-jr ) = y^sinfc^ sing y (g: k - £ k+Q Jc t (q)c(q + Q y )c f (k)c(k + Q 3 

k,q 



dt 



dt 



- (x <r+ y) 



(2.30) 



Now we can factorize the 4- Fermi term using (c^(k)c(k)) = n(k): 

= ^ sin k x sin k y c j (k) c(k + Q) (e k+Qy - e k +q)(l - n(k + Q v )) 



dt 



dt 



- (e k - e k+Qx )n(k+ Q x ) - (x y) 



sin k x sin fc y c^(k)c(k + Q) 



+ 2(e k+Qx n(k + Q x ) - e k+Qy n(k + Q y )) 
- (e k + £ k +Q)(^(k + Qx) - n(k + Q y )) 



(2.31) 



The r.h.s. is indeed ~ Q. 



III. FLUCTUATING NEEL AND VBS ORDERS 

Given the connection between the skyrmion number of the Neel order and the VBS order 
derived in Section [TTJ it is natural to wonder whether the two order parameters can be 
treated at a more equal footing. In Section [IT] we investigated the fermion correlations in the 
background of a spatially varying Neel order, and so this suggests a natural generalization 
in which we allow for a background spacetime dependence of both the Neel and VBS orders. 
This section will present the needed generalization. The result here will be an alternative 



derivation of the arguments of Section II B : the skyrmion number of the Neel order and the 
angular variable, 0, of VBS order V ~ e 1 ^ are quantum-mechanically conjugate variables. 
We start from a Neel state with the order parameter m a = m(n x , n y , n z ) ^ 0. When the 
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system approaches the Neel-VBS transition, fluctuating VBS order becomes important and 
needs to be taken into account. The starting point of our analysis is the electron Hamiltonian 
H with both m a and V: 

H(V X , V y , n x , n y ) = [H b + mH?] + [V X H^ + V y E v y + m(n*H% + n v H?)) (3.1) 

where iff, describes the electron band structure in the absence of Neel or VBS order as in 



Eq. (2.1); the fermion bilinear operator (H^,Hy) is dimerized electron hopping in x and 



y directions; (H^ , Hy , Hf) is staggered electron spin density in x,y,z directions in spin 
space; (V x ,V y ) describes the fluctuating VBS order; (n x ,n y ) describes the Goldstone mode 
of the Neel order. We now integrate out the fermions and derive the effective action S for 
the slowly varying bosonic fields A^(x, y, r) = (V x , V y , n x , n y ). 

Treating the second term in if as a perturbation, we find couplings between Neel and 
VBS order starts at fourth order in a one-loop expansion: 

r 3 

5l = E / U d ^ K ^L-^(PM u (P2)A x (p 3 )A s (-p 1 -p 2 -p 3 ), (3.2) 

[Note: we drop all numerical prefactors in this subsection.] Here fj,,v,\,5 = 1, ...4 labels 
the components of the perturbation field and the vertex; p = (p°,p x ,p y ) is the external 
momenta of A^. We now expand the function K in powers of p and collect terms involving 
the product vlv^vi with a, {3, 7 = 0,x,y: 

k^L = k;^- v ^vI + ... (3.3) 

This corresponds to a derivative expansion in real spacetime: 

Si= Yl K Zi & J dxd V dr {A s d a A»d[ ) A v d r A x ) (3.4) 



The action (3.4) resembles the Chern-Simons theory in 6+1 dimensions. A difference is 
that the space-time indices a, f3, 7 and the internal indices /1, u, 5, X do not mix with each 
other. Qi et al^ recently proposed that Si can be simply obtained from the Chern-Simons 
term by dimensional reduction to 2+1 dimensions. The procedure is to throw away all 
components in the Chern-Simons term, which involve spatial derivatives in the internal 
dimension. We shall show by calculating K^ x,s explicitly that this dimensional reduction 
approach does not apply in the present situation. 

Among the terms in S\, we are particularly interested in a topological term 



S top = ^2 / iK a p J dxdydrj^j^, 



(3.5) 
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where j£ is the skyrmion current in the Neel state: 

3a = e a /3~ / eabcn a d[ 3 n b d~ / n c , (3.6) 

and jj is the VBS current: 

jj = V x dpV y - V y d p V x . (3.7) 

It follows from symmetry analysis that on the square and honeycomb lattice, The matrix 
K a p is diagonal. Because of four- and six-fold rotational symmetry, K xx = K yy . S top then 
becomes 



S top = i J dxdydr (Kj?jY + .j]' + iCtfg) 



(3.8) 



Comparing (3.8) and (3.4), we can express K in terms of the tensor components 



K oc [K Qxy ' + Permutations of (2,0), (3,x), and (4,?/)] 

— [Kq^' 1 + Permutations of (2,0), (4,x), and (3,y)] 

— [K^y' 2 + Permutations of (1,0), (3,x), and (4,?/)] 

+ [K^y 2 + Permutations of (1,0), (4,x), and (3,y)] (3.9) 



We now calculate K for the square lattice. The Hamiltonian Hi, is specified in Eq. (2.1), 
and we choose only nearest neighbor hopping t. For the coupling to the order parameters, 
we choose 



#f = 5>t( ri Kc( ri ) - JVfoKcfc) 

Hj = ^(-l) i "[ C t(r i )c(r i + e i9 ) + c.c.], /3 = x,y (3.10) 

where we have divided the square lattice into two sublattices A and B defined by (— iy x+l y = 
±1. The Neel order carries crystal momentum (vr,7r). The VBS order in the x- and y- 
directions carries crystal momentum (jr, 0) and (0, tt) respectively, with the corresponding 
dimerization pattern shown in Fig. [3} Note that the dimerization pattern "rotates" around 
a site as the phase of V x + iV y advances by 2ir. It is straightforward to check that the term 



Stop in Eq.(3.4) satisfies square lattice symmetry. 

The lattice periodicity is doubled in both the x and y direction for m ^ and V x , V y ^ 0. 
We choose the 4 sites in a plaquette as the new unit cell. The Bloch Hamiltonian H(k x , k y ) 
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FIG. 3: VBS order on the square lattice, 
is obtained by Fourier transform: 



H(k x , ky 



/ 


m a a a 









\ 




t* 
l y 


—m a a a 











t* 





—m a a a 


y 




\ 





t* 


t* 


m a a a 


J 



-t cos ky + iVy sin k y , 
-t cos k x + iV x sin k x . 



(3.11) 



K = 



Here k x ,k y G [— 7r/2,7r/2] is crystal momentum in the folded Brilluoin zone. We used 
Mathematica to evaluate K and found 

" dk dk dk ra 3 £ 2 sin 2 k x sin 2 ^[t 4 (cos 2 k x - cos 2 k y ) 2 - (fcg + m 2 ) 2 ] 

x y [t 2 (cos k x - cos k y ) 2 + k% + m 2 } 3 [t 2 (cos k x + cos k y ) 2 + k% + m 2 } 3 ' 

The integration over ko can be done analytically using Mathematica. The resulting integrand 
K,(k x , k y ) is a complicated function of k x and fc r Instead of showing its explicit form, we 
plot )C(k x , k y ) over the Brillouin zone k x , k y £ [0, ir] in Fig. [1] Note that the integrand is 
peaked at the "hot spot" Q = (7r/2,7r/2). This is not surprising because both the Neel and 
VBS orders have strong nesting at Q. 
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FIG. 4: Plot of JC over a quarter of the Brillouin zone for t = 1, m = 0.25. 



The other coefficient K' in Si can be obtained similarly and is given by: 



K' 



dkodk x dh, 



mH 2 sin 2 k x sin 2 k y [t 4 (3 cos 2 k x + cos 2 k y )(cos 2 k x + 3 cos 2 k y ) — 3 (A; 2 . + 



\t 2 (cosk x — cos k y ) 2 + k 2 + m 2 ] 3 [t 2 (cos k x + cos k y ) 2 + k 2 , + m 2 } 3 
Comparing K and K', we found that in general K ^ K'. This means that different terms 



in the effective action (3.4) have different coefficients, so that S cannot be obtained by 



dimensional reduction from a Chern-Simons term in 6+1 dimensions, which has a single 
coefficient. 



A. Quantum disordered Neel phase 

We can now use the results of this section to present an alternative version of the argument 



in Section II B that the quantum disordered Neel phase has VBS order. The argument here 
will be closer in spirit to the duality mapping discussed in Ref. El 

We will limit our discussion to a quantum-disordered Neel phase where the monopole 
density is very dilute. Thus, we assume that over a significant intermediate length scale 
there is an effective description in terms of a theory in which the total Skyrmion number is 
conserved. As discussed in Section [TT] and Appendix [Bj we can represent the fluctuations in 
the local Skyrmion density by a low energy U(l) photon field A a : by Eq. (2.7), the gauge 
flux in this field, e a prydpA 7 , is a quarter the Skyrmion current j£ in Eq. (3.6). We can write 
an effective action of the photons as 



C, 



eff 



1 

2^2 



(3.13) 



Here the second term represents the topological term in Eq. (3.8). For simplicity, we have 
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assumed K' = K. Different values of K and K' will not affect our conclusion below. Also 
note that by the discussion at the end of Section II B Q is conjugate to 



Now let us perform the standard duality transformation of 2+1 dimensional 
electrodynamics^^ on £ e g. The first step corresponds to decoupling the Maxwell term 
by a Hubbard-Stratonovich field, J a , to obtain 



£eff — TT^Ja + 7T~ Ja^aB-ydpA^ + 4AK € a /3~,dpA~ . 

Now, we integrate over A a , and this yields the constraint 



Ja = d a (f> — 8nK j 



v 



(3.14) 



(3.15) 



where is the scalar field which is dual to the photon. We have judiciously chosen factors of 
(2?r) above to ensure a normalization so that e 1 ^ is the monopole operator. Finally, inserting 
Eq. ( |3.15 ) into (3.14) we obtain 



C eS = ^(d a <j>-8<KKjV)' 



(3.16) 



The effective Lagrangian for the photon phase in Eq. (3.16) allows to conclude that the 
long-range correlations of d a cf> have the same form as those of j^. In other words, we have 
the operator correspondence d a cj) ~ j^. In terms of the complex VBS order parameter 
V = V x + W y we can therefore write for the monopole operator e 1 ^ ~ V u , where v in general 
appears to be an irrational number. In the special cases where the value of K was quantized 
by projection to an integer number of electrons per site^^, v was found to be an integer; 
this is a possible value of v here, although our present methods don't allow us to see why any 
particular integer would be preferred. The uncertainty in the value of v here is analogous 
to the arbitrariness in the overall normalization of Q in Section |II A 



In any case, as long as v is not an even integer, the correspondence between the monopole 
operator e 1 ^ and V u implies that V has long-range correlations in the monopole-free region. 
At even longer scales, once the monpoles condense, the phase of V is locked along one of 
the lattice directions 89 . 



IV. HONEYCOMB LATTICE 

This section will apply the methods developed in Section [IT] and Appendix [B] to the 
honeycomb lattice. As is well known, this lattice has an electronic dispersion with a Dirac 
form at low energies. We will adapt our methods to the Dirac fermions, and find that many 
results can be computed rapidly in closed form. 

The honeycomb lattice has 2 sublattices, and we label the fermions on two sublattices as 
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ca and c B . To begin, we only include Neel order explicitly. Then the analog of Eq. (2.2) is 



H = -t {c Ai c Bj + c Bj c A ?j + m Y c ] Ai n a {Yi)a a c Ai - m ^ c t m n a (r i )a a c Bi . (4.1) 



(«■> 



We restrict to the case with constant Neel order n a , transform to momentum space, and 
introduce Pauli matrices r a in sublattice space, and obtain 

/d 2 k \ ( \ 
— c f (k) -t(cos(k ■ ei) + cos(k • e 2 ) + cos(k ■ e 3 )jr x 

+ t ( sin(k • e x ) + sin(k ■ e 2 ) + sin(k ■ e 3 )) r y + mr z n a a a ] c(k) (4.2) 



(4.3) 



where we have introduced the unit length vectors 

e 1 = (l,0) , e 2 = (-l/2,V3/2) , e 3 = (-1/2, -^3/2). 



We also note that we take the origin of co-ordinates of the honeycomb lattice at the center 
of an empty hexagon, so the A sublattice sites closest to the origin are at ei, e 2 , and e3, 
while the B sublattice sites closet to the origin are at — e 1; — e 2 , and — e 3 . 

The low energy electronic excitations reside in the vicinity of the wavevectors ±Q 1; where 
Qi = (47r/9)(e 2 — e 3 ). So we take the continuum limit in terms of the 8-component field C 
defined by 



C 



Al 



ca(Qi 



c 



Bl 



c B (Q 



1 J 



C 



A2 



c A (-Q!) , C B2 = c B (-Q 1 ). (4.4) 



In terms of C, we obtain from Eq. (4.2) 

72 



H = J ^C f (k) (vT y k x + VT x p z k y + mT z n a a a ^C{k), 



(4.5) 



where v = 3t/2; below we set v — 1. We have also introduced Pauli matrices p a which act 
in the 1, 2 valley space. This is the final form of H: it makes the Dirac structure evident, 
and will also be the most convenient for our computations. 

It is also convenient to list the effects of various symmetry operations on C . Under 



reflections, I y , which sends x -H- — x 



'■ Cai — > C B \ 



Al 



Cai — >• C 



B2 



Cro — > C 



A2 



(4.6) 



Similarly 



T-x '■ Cai — > C 



'A2 



C m -> C 



B2 



Ca2 — > C 



Al 



Cm — > C 



Bl 



(4.7) 
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FIG. 5: VBS order on the honeycomb lattice. 



Rotations by 60 degrees, R, lead to 

R : Ca\ — > u 2 Cb2 , Cbi uC A2 , Ca2 uCbi , Cb2 u 2 Cai (4.8) 
Translation by the unit cell vector e 2 — e3, T y : 

T y : Cai — > uj 2 Cai , Cbi u 2 Cbi , Ca2 uCa2 , Cb2 wC B2 (4.9) 
Time reversal t — > —t: 

T: C A1 ^ia y C A2 , C A2 ^ia y C Al , C B1 ^ia y C B2 , C B2 ^ia y C m (4.10) 

Notice that time reversal transformation also involves a complex conjugation transformation. 

From these transformations, we can construct the fermion bilinear associated with the 
kekule VBS pattern shown in Fig [5} In terms of the continuum field C, the VBS order 
parameter is 

V = C ] t x ( P x + ip y )C (4.11) 

We can verify this is the VBS order with the kekule pattern of Fig. [5] by its symmetry 
transfer mat ions 



T 

■hi 


V ->■ V 


T 


V -> V* 


R 


v^v* 


T 


V u 2 V 


r 


V V*. 



(4.12) 



A. 6D method 



In the present situation with a Dirac fermion spectrum, the dimensional reduction 
methocP^l from 6D does apply, and be used to compute the coupling between the fluc- 



tuating Neel and VBS orders. From Eqs. (4.5) and (4.11), we can write down the explicit 
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form of the Hamiltonian in the 8x8 space of Dirac fermions: 



H(k, ki, k 2 , h, k 4 ) = 

r y k x + T x p z k y + T x p x h + r x p y k 2 + r z a x h + r z a y h + mr z a z n z . (4.13) 

Here k = (k x ,k y ), and the 'extra-dimensional' momenta £4,2,3,4 are related to the order 
parameters: k\ = V x , k 2 = V y , k% = n x , and k^ = n y . Now note that the matrices in all the 



terms in Eq. (4.13) anti-commute with each other. So this has the natural interpretation 
as 6D Dirac Hamiltonian, where the last term proportional to m has the interpretation of a 
Dirac fermion mass. We can now proceed as in Ref. UHl and derive the WZW term for the 
order parameters. 



B. U(l) gauge theory 



Next we turn to the analog of the analysis in Section [HA] for the square lattice. However, 



rather than working with the spatially varying Neel order as in Eq. (2.2 ), we will use the U(l) 



gauge field formulation of Eq. (2.5) which was applied to the square lattice in Appendix [B 



We begin with the ip fermion Hamiltonian in Eq. (2.5) and take its continuum limit as in 



Eq. (4.5). For this, we make the analog of the transformation in Eq. (4.4) from the lattice 



ifj fermions to continuum \1/ fermions. In this manner, we obtain the continuum U(l) gauge 
theory 



C = ^{d T - ia z A T - ir y (d x - ia z A x ) - ir x p z (d y - ia z A y ) + mr z a z ^ 



(4.14) 



Now we obtain the result which is the analog of Eq. (BIO) by applying the Kubo formula 



to determine the response in an arbitrary fermion bilinear ^/''T^/ due to an arbitrary slowly 
varying A a . This involves evaluating a diagram with one fermion loop, and the long wave- 
length result is 

<^r^) = ^-^d x A T -d T A x )TT[TT x ] + (d T A y ~d y A T )Ti[Tp z r y ] 

+ (^.-^)Tr[r P 2 ]). (4.15) 



Now we see that the choices T = t x , p z r y , and p z lead to non-zero fermion bilinears induced 
by the A a gauge flux. Note that this result was obtained with much greater ease in the 



continuum Dirac theory than for Eq. (BIO) 



Let us restate the result in Eq. (4.15) in different terms. We add to C in Eq. (4.14) a 



source term j^: 



(4.16) 
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Then the implication of Eq. (4.16) is that if we integrate out the ¥ fermions, the effective 
action for and the gauge field A a has a mutual Chern-Simons term: 



£ 



off 



1 2 * V 



(4.17) 



The similarity to Eq. (3.14) should now be evident. We can now proceed with the duality 



of electrodynamics to obtain the analog of Eq. (3.16), which is 

3m 



£ 



eff 



•>7T 



(<9 Q - jl) 



(4.18) 



where again e 1 ^ is the monopole operator. As argued below Eq. (3.16) any operator with 



the same quantum numbers as e 1 ^ has long-range order in the 'quantum-disordered' phase. 
Here we present a different route to identifying candidates for the competing order. First, 



we notice that the theory in Eq. (|4.16|) actually enjoys a gauge invariance under which 

(4.19) 



.f>- 



¥ -+ exp ( ) ¥ , f a ^jl-d a 9 



where 9 is a field with an arbitrary spacetime dependence. (Note that this gauge invariance 
is completely different from that associated with the A a gauge field, under which ¥ — > 



exp((i/2)a z 9')ty.) Now we observe that this gauge invariance extends also to Eq. (4.19) 
under which 



(4.20) 



We will use Eq. (4.20) as the key relation needed for any competing order associated with 
the monopole operator e 1 ^. 



Equivalently, we can use Eq. (4.19), and restate the requirement of Eq. (4.20) as the 
commutation relation 

[Q,e^ {x) ] = e^ (x) , (4.21) 



where 



Q=- I d 2 r $V$. 



(4.22) 



This makes a very explicit connection to Section II B and Eq. (2.28 ). Note that here the over 



all normalization of Q is specified, and does not suffer from the arbitrariness we encountered 
in Sections III Al and IIHAI 



Now we can easily check that the VBS order parameter in Eq. (4.11) obeys the commu- 
tation relation 

[Q,F(x)] = F(x), (4.23) 
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and so we conclude that e 1 ^ ~ V, and that VBS order can appear in the quantum-disordered 
Neel phase. 



C. Other competing orders 

In addition to the VBS order parameter V, it is now easy to see that there are other 
order parameters which are canonically conjugate to Q. For instance, the following three 



complex order parameters all satisfy Eq. |4.23 

Vi ~ &(p x + ip y )ty, V 2 ~ &T z (p x + ip y )V, V 3 ~ ¥r y (p x 
Under discrete symmetries, these order parameters transform as 



ip v )^). 



(4.24) 



T 

■hi 


Vi^V x , V 2 ^-V 2 , V 3 ^-V 3 , 


(4.25) 


T x 


v, v; 


(4.26) 


R 


Re[V]] + ilm[V 2 ] -> u; 2 (Re[^] + ilmfV^]), 


(4.27) 




Re[V 2 ] + zIm[Vi] ->• -u 2 (Re[V 2 ] + ilm^]), 


(4.28) 




v 3 -> -v* 


(4.29) 


T 
± y 




(4.30) 


T 


Vi v* , v 2 ^v* , v 3 ^-v;. 


(4.31) 



According to these transformation laws, we can identify that V\ is a charge density wave 
(CDW) with wave vector 2Qi, V 2 is the A — B sublattice staggered CDW, and V 3 is a charge 
current density wave. 



However, notice that the matrices in V in Eq. (4.11 ) anticommute with all the matrices in 



H in Eq. (4.5); therefore the VBS state has the lowest fermionic mean field energy, because 
the fermion \1/ will acquire a Dirac mass gap m ~ ^m 2 + |l^| 2 . Compared with the VBS 
order parameter V, the other three order parameters have higher mean field fermion 
energy, hence are less favorable in energy. 



D. Superconductor order parameters 



In addition to the VBS order parameter, the Neel order can also have strong competition 



with superconductor, as long as the SC order parameters satisfy Eq. 4.21. In this section we 



will focus on spin singlet pairings. Using the quantum number Q in Eq. |4.22 and criterion 



Eq. |4.21[ it is straightforward to show that the following six groups of SC order parameters 
are candidate competing orders of the Neel order: 



Group 1 : (A x , A 2 ) ~ (Re[VHa y % Imbl^V*]), 
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Al ~ ^ C , A,Q 1+ fc^Cr 2/ CA,Q 1 -fc + CB j Q 1+ fci<7 S/ C.B,Q 1 -fc 

fc 

+C.A -Qi+fcic^CU ,-Qi-fc + C'fi -Q 1+k i<7 y C B ,-Qi-fc + # c. 
A 2 ~ ^«CU,Qi+fc«C^CU,Qi-fc + iCs.Qi+fc^Cs.Qi-fc 

fc 

-iC A -Q 1+k ia y C A -Qj-fc - iC B -Q 1+k ia v C B -Qi-fc + #-c 



Group 2 : (Ai, A 2 ) ~ (Im[$^^], Reftf'icrV*]), 

Ai ~ ^ »CU,Qi+fc^cr y CU,Q 1 -fc + iC , B) Q 1+fc i<7 s/ Cs j Q 1 _fc 

+*Ca -Q 1+k ia y C A -Q^fc + iC B - Ql+k ia y C B -Q 1 - k + if.c. 

A 2 ~ ^ Cyt,Qi+fc»cr i/ C > A,Qi-fc + C , B i Q 1+ fcicr 2/ C S) Q 1 _fc 
fc 

— CU-Qi+fcic^CU-Qi-fc — Cs-Qi+fc^^Cs-Qi-fc + -^- c - 



Group 3 : (Ai, A 2 ) ~ (Re^VW*], Im[*V z i(rV*]) 

Al ~ ^ C > A,Q 1 +fci(7 i/ C f A,Q 1 -fc - Cs^+ki^CB^-k 

k 

+C A -Q 1+k ia y C A -Qi-fc - Cb -Qi+fc^^Cfi ,-Qi-fc + -^- c - 

A 2 ~ ^ iC A ,c tl+k ia y C A Q. L - k - iC B ,Q 1 +kicr v CB,Q 1 -k 
k 

-iC A -Q 1+k ia y C A -Qj-fc + iC B _Q 1+ fcicr y C B _ Ql _fc + -ff.c. 



Group 4 : (Ai, A 2 ) ~ (Im[^W#], ReftfViaV*]) 

Ai ~ ^ »CU,Qi+fc^ y CU,Qi-fc - ^C^Qi+fc^C^Qi-fc 

fc 

+iCU -Qi+fcio-^CU ,-Qi-fc - iC B -Q l+k ia y C B -Qi-fc + #.c. 

A 2 ~ ^ CA,Qi+fc^ y CA,Qi-fc — C^Qi+fc^C^Qi-fc 

fc 
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— Ca -Q 1+ kiCT y CA ,-Qi-fc + Cb -Qi+fe^^CB-Qi-fe + ^ -C 

Group 5 : (Ai, A 2 ) ~ (ReftfVWtf], Imft* 7* ia v ffV]) 

Ai ~ ^ CA.Qi+fc^Cg.Qi-fc + Cs^+ki^CA^-k 
k 

+C A -Q 1 +kicr y CB-Q 1 -k + Cb -Qi+fcicr^CU ,-Qi-fc + if.c. 

A 2 ~ ^ iCA^+ki^CB^-k + iCB&j+kio-vCA&^k 
k 

-iCA-Qi+ki^CB-Qi-k - iC B -Q 1+k icr y CA -Qi-fc + #c. 

Group 6 : (A 1; A 2 ) ~ (ImftfVWtf], Re[^Wp 2 ^]) 

Ai ~ ^ iCA^+ki^CB^-k + iCB,Q 1 +kicr v CA,Q 1 -k 
k 

+iC A -Q 1 +ki(7 v CB-ci 1 -k + «C B _Q 1+k ia y C A ,-ch-fc + #c. 

A 2 ~ ^ C^Qj+fc^C^Ch-fc + CB,Q 1 +A,i0" i 'CA,Q 1 -fe 
fc 

-Ca -Q 1+k ia y C B -Qi-fe - C B -Q 1+k i<j y C A -Qi-fc + #c. (4.32) 

All of these SC order parameters carry nonzero lattice momentum 2Qi, and none of them 
gaps out the Dirac points. Nevertheless, these SC orders are most likely to be adjacent to 
the Neel order on the phase diagram. 

V. NAMBU QUASI-PARTICLES OF d-WAVE SUPERCONDUCTOR 

In this section we will apply the above methods to analyze the d— wave superconductor 
and its descendants. As in previous section we will examine the nature of the "quantum 
disordered" phase after loss of antiferromagnetic order. However, we will not consider the 
case of commensurate antiferromagnetic ordering at wavevector Q = (jr, ir), because it 
requires computations we have not explored here. Rather, we will limit ourselves to the 
technically easier case of nested spin density wave order, with a wavevector precisely equal 
to the separation between two of the nodal points of the fermionic excitations of the <i-wave 
superconductor. The non-nested case is of experimental importance, but we will not consider 
it here. 
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The nodal quasi-particles of the d-wave superconductor are described by the following 
Dirac fermion Lagrangian: 

= y\(d T -t^(d x + d y )T Z -t^=(-d x + d y )T^ 1 

+ *l(er ~ ^(- d * + d vV - l ^ d * + dy)r X )^ 2 . (5.1) 

^1 = (fi, i^flY, ^2 = (f2,icr y flY- fi, /2, /3 and f± are quasiparticles at nodal points 
{Q,Q), {—Q,Q), {~Q,—Q) and (Q,—Q) respectively. Notice that Q is in general incom- 
mensurate, vp and i>a are different from each other. 

We assume the system has the symmetry of the square lattice. Under square lattice 
discrete symmetry group, the quasi-particle and transform as: 

T x , x ->■ x + 1, *i ->• e lQ V u ^ 2 ->• e- lQ ^ 2 ; 
T y , y -> y + 1, *x -)• e^, ^ 2 -)• e lQ ^ 2 ; 

I x , 2/ ->■ -y, ^1 ->■ cr y r y ^l, ^ 2 ->■ crV^l; 

2^, x ->■ y, y ->■ x, *x ->■ zt**i, ^ 2 -)• txVttJ; 

X, +2/ , x->-y, -x, ^1 ->• aV^l, ^ 2 ->■ ir*tf 2 ; 

T, t ->■ -i, ->■ ir y ^\, ^ 2 ->■ ir w *5. (5.2) 

Notice that transformations T x _ y and X x+y are combined with a U(l) transformation on the 
superconductor order parameter: Ae l9 — > Ae te+tn . 

Now let us consider the spin density wave states that gap out the nodal quasi-particles i.e. 
the SDW with wave vector (2Q,2Q), and (2Q, —2Q), which can be written as m>\r y a y a^> 1 
and i^! t 2 T y a y 2 respectively. In contrast, the SDW at (2Q, 0) and (0, 2Q) will not gap out 
the nodes, and they will be ignored hereafter. It is convenient to introduce the Majorana 
fermion Xa as \& = xa + ^Xb, and there are in total four different choices of SDW that can 
gap out the nodal points: 

$! = Re[^WW] ~ (x*irWp 2 x, x t r y (T y a y p x X , X*irW 'p z X ), 

$2 = Be[M t T v a y atJi z *] ~ (x*irWpVx, x^WpVx, x'^WpVx), 
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$ 3 = Im[^VW#] ~ (x t ir y a y a x p x X , xVctVVx, X^r y (J y a z p x x), 

$ 4 = lm[i¥r y a y afi z ^\ ~ ( X t ir y a y a x p x p J z X , x'rWpVx, X*zrWpVx)- (5.3) 

The Pauli matrices p a mix and \l/2, while the Pauli matrices p a mix ^ and X b- In the 
Majorana Fermion basis, the SU(2) spin operators are represented by the total antisymmetric 
matrices 

S= (a x p y , a y , a z p y ). (5.4) 

— * — # 

We can check that all four vectors $ a (a = 1 • • • 4) transform as vectors under S. 

Mow we hope to consider the slowly varying SDW by introducing the SU(2) gauge field 
Y^i A^S 1 , which will be Higgsed to U(l) gauge field with a background nonzero expectation 
value of $ a : 

L x = x\((d T - iA l S l ) - iv F (d x - iA l x S l )r z - iv A {dy - iA l Y S l )r x ) X i 
+ xt({d T - *A l S l ) - iv F (d Y - iA l Y S l )r z - iv A (d x - iA l x S l )r x ) X 2 

+ KxXx- (5.5) 

Now we have redefined the coordinate (x + y)/y/2 — > X, (—x + y)/y2 — > Y. The order 
parameter $^ has Higgsed the SU(2) gauge field down to U(l) gauge field A 1 . The matrix 



T a can be found in Eq. 5.3 



Now the flux quantum number can be calculated using the same techniques developed in 
the previous sections. We summarize our results in the following: 



Group 1 


($i> * o, 


gauf 


?e flux carries Q ~ 


X t cr y p z p x X ; 


Group 2 


(<H> * o, 


gauf 


ie flux carries Q ~ 


X t v v P x X, 


Group 3 


<$ 3 > ^ 0, 


gauf 


ie flux carries Q ~ 


X l (r y p z p z x, 


Group 4 




gauj 


$e flux carries Q ~ 


X^ y p z X- 



(5.6) 

The quantum number Q carried by the flux is obviously SU(2) gauge invariant. 

The flux condensate will again lead to orders that break certain lattice symmetry in 



Eq. 5.2 The condensate order parameter V has to satisfy Eq. 2.28 Within all these order 



parameters that satisfy Eq. 2.28, we choose the order parameters that have the lowest nodal 



quasi-particle mean field energy, i.e. the order parameters that anticommute with T l a . We 
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list our results in the following equation, and for each group of SDW in Eq. 5J3 we introduce 
a five component vector E\ a) with $ a ~ (Sj o) ,S 2 a) , Sj a) ), V a ~ Hj a) + iSj ft) : 

Group 1 : H/^ ' ' 3 = $i, 



Group 2 : EL, ' 2 ' 3 = $ 2 , 



Group 3 : HLj ' 2 ' 3 = $ 3 , 



Group 4 : SjT, 1 ' 2 ' 3 = <£ 4 , 



(4) 

1 

-(4) " -(2) " ^ 



1 



H* 4) = E 5 {1) = — X \t z - r*VWx ~ Imf^r* - r^^^]. (5.7) 

With the formalism developed in Ref.^, we can also show that there is a 0(5) WZW term 
for each group of 0(5) vector Ely Both the WZW term and the gauge flux calculations 
imply that the SDW $| s and order parameters V a ~ Ei^ + zHjK are competing with each 
other, and after suppressing the SDW the system enters the order with nonzero (V a ) 
directly. 

Now we want to identify the physical meanings of H?n and SL. Clearly, Six and H^U are 
both density waves of physical quantities, with wave vectors (2Q, 0) and (0, 2Q) respectively. 
Under lattice symmetry Eq. 



5.2 



and EjK transforms as: 



•7- . r- >4 , '4 "5 . "5 "4 , ^4 "5 . ^5 
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• T 

->-x 


—4 

"(1) 


— > 


"4 

"(1)' 


"5 


"5 

"(1) 


x 

-^x—y 


"(1) 




^5 

"(2)' 


"(1) ° 


^4 

"(2) 


1-x+y 


—4 
-(1) 


•H- 


■^5 
"(2)' 




^4 

"(2) 


T 


S (a) 


— >■ 


S (a), 


i = 4, 5 





"(2) -> "(2), "(2) -> --(2), 



(5.8) 

According to these transformations, we can make the following identifications: 

5j 2) + iSj 1} = Sj 4) + z'Hj 3) = VBS or CDW with wave vector (2Q, 0); 

Sj? 1} + z'E^ 2) = S^ 4) + z'Sj? 3) = VBS or CDW with wave vector (0, 2Q). (5.9) 

These analysis suggests that the SDW at wave vectors (2Q, 2Q) and (2Q, —2Q) is competing 
with CDW/VBS order parameters at (2Q,0) and (0,2(5), an d the suppression of the SDW 
leads to the emerging of CDW/VBS order parameters. 



VI. CONCLUSIONS 

This paper has addressed a problem of long-standing interest in the study of correlated 
electron systems in two spatial dimensions. Many such systems have insulating, metallic, 
or superconducting ground states with long-range antiferromagnetic order. By tuning the 
electron concentration, pressure, or the values of exchange constants in model systems, it is 
possible to drive a quantum phase transition to a phase where the antiferromagnetic order 
is lost. We are interested in the nature of the "quantum-disordered" phase so obtained. 

For certain insulating square or honeycomb lattice models, the essential features were 
understood some time agcP® the lattice spins endow point spacetime defects in the Neel 
order ('hedgehogs') with geometric (or Berry) phases, which lead to valence bond solid (VBS) 
order in the quantum- disordered phase. Here we have presented a more general version of 
this argument, in principle applicable to arbitrary insulating, metallic, or superconducting 
electronic systems in two dimensions, with general band structures. The key step was 
to associate the geometric phases with bands of one electron states in the background of 
local antiferromagnetic order. The antiferromagnetic order was then allowed to have a 
spacetime variation in orientation (but not in magnitude) so that there was no long-range 
antiferromagnetic order, thus accessing the quantum-disordered phase. We found that the 
skyrmion density in this local antiferromagnetic order induced a response in an electronic 



bilinear conjugate to the competing order: this is contained in our key result in Eq. (2.22). 



Our main application of these results was to cases in which the electronic band structure 
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was fully gapped in the phase with antiferromagnetic order: we considered square lattice 
insulators in Section [TTJ honeycomb lattice insulators in Section IV, and (i-wave supercon- 
ductors with spin density wave order nesting the nodal points in Section V We obtained 
VBS order in many cases, but also found a number of other possible orderings. 

However, in principle, the result Eq. (2.22) applies also in cases where the antiferromag- 
netic order does not fully gap the electron bands e.g. when there are hole and/or electron 
pockets. Such a situation is clearly of importance for the underdoped cuprate superconduc- 
tors. The result in Eq. (2.22) contains a singular dependence on k at the Fermi surfaces 
of such band structures, and this is likely of importance in the quantum- disordered phase. 
Alternatively, expressions for the coupling K in Section III would acquire long-ranged cor- 
rections due to Fermi surface singularities. We leave the elucidation of such effects to future 
work. However, if we ignore such effects, the arguments of Section HI A would apply also 
to this metallic case, with a variable exponent v relating the monopole operator to the VBS 
order. The net result is that any ordering associated with an integer power of V is possible. 
Interestingly the same conclusion was reached in an earlier study®! Q f quantum disordered 
Neel states in a compressible background using a toy model of bosons. 
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Appendix A: Rotor theory of Hubbard model 



This appendix will show how the decomposition in Eq. (2.3) can be used to write an 



exact path integral representation of an arbitrary Hubbard-like model. The z a becomes 
co-ordinates of an 0(4) rotor in this path integral, and so do not directly contribute to the 
geometric phases of interest in this paper. This is to be contrasted from the alternative 
Schwinger boson formulation, where the canonical nature of the Schwinger bosons ensures 
that they carry the entire geometric phase at half-filling^. 
We consider a Hubbard model on a general lattice 



H = H + H 1 



(Al) 



where H Q has the single site terms 



U n it 



(A2) 
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and Hi is the hopping term 



(A3) 



As in Eq. (2.3), we transform the electron to a rotating reference frame expressed in terms 
of the spinless fermions c p and the complex unit spinor z a . Here, it is useful to write z a in 
real and imaginary parts: 



Z\ = 00 + i<f>l , H = 02 + 203- 



The inner product of two complex spinors is 



r a z a = 0(1 - p y )<p 



(A4) 



(A5) 



We will use a a for Pauli matrices in the f, I space, and p a for Pauli matrices in the 
real/imaginary space. The global spin rotation 



(l + i9 a a a )z 



acts on via 



0^ (l + i6 a S a ) 

where S a are the antisymmetric Hermitian matrices 



S x = -a x p y , S y = a y , S z = -a z p y . 



Combining (A5) and (A6) we have 



z *a a z = 0(1 - p y )S a 4> = -4>p y S a 4>. 
The SU(2) gauge rotatiorP acts on ip as 

V> ->■ (1 + i6 a a a )4> 

where a a are Pauli matrices in the ± space. This gauge rotation acts on z as 



-)• (1 + OT" 



£m7 rim 



(A6) 



(A7) 



(A8) 



(A9) 



(A10) 



(All) 



where the indices £,m — 1 ... 4 and T a are the antisymmetric Hermitian matrices 



T x = a y p x , T y = -a y p z , T z = p y . 



(A12) 
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The physical states on a single site, which are eigenstates of H 0} are 



4|0) ^ (zWl-eapz^-jlO) 
|0> |0) 
cjc{|0) <-> V'+V'llO) 



(A13) 



These 4 states have energies — // — 11/4, —ji — U/A, Z7/4, and Z7/4 — 2/i. 

Following Hermele®', let us write these states in a different manner, using the energy 
levels of a 0(4) quantum rotor. All of the following will work on a single site, and so we will 
drop the site index. We will equate the states of <ft to that of a quantum particle moving on 
S 3 with co-ordinate (j). On this space, we introduce the angular momentum operators 



d 



d 



(A14) 



In the fermion sector we have the usual angular momentum 



Then all the states in Eq. (A13) satisfy 



T a + C a = 0. 



Now consider the following Hamiltonian for the rotor and the fermions 
Ho = K ± S a2 + K 2 T a2 + K 3 ^ p + K^U_^, + + K 5 



(A15) 



(A16) 



(A17) 



For appropriate ranges of the Ki couplings, the low-lying states of this Hamiltonian which 



obey Eq. (A16) map onto the states of the Hq. The zero rotor-angular momentum states 



must have or 2 fermions, and these map onto the lower two states in Eq. (A13), yielding 



K 5 

2K 3 + K A + K 5 



U/A 

U/A - 2fi 



(A18) 



There are 4 rotor states with angular momentum 1 and wavefunction ~ Because of 



the constraint in Eq. (A16), these states must be paired with states with fermion number 1. 



There are 2 such states, leading to a total of 8 states. However, the conditions in Eq. (A16) 



eliminate 6 of these states (there are 3 constraints for each fermion polarization), and so 
only 2 states remain, as in the Hubbard model. The energy of these states yields 



3#i + 3K 2 + K 3 + K 5 = -fi- U/A. 



(A19) 
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The Kj constants are over-determined, and in an exact treatment of the constraint in 



Eq. (A16), the precise choice will not matter. Of course, in mean-field theory, different 



choices will lead to somewhat different results. 

Now, following HermekP^, we can write Eq. (A17) as a path integral over 4>t(j) and ip p {r) 
and obtain the Lagrangian 



1 



4(^i + K 2 ) 



(A20) 



where A%. is the time-component of a SU(2) gauge field which imposes the constraint (A16), 
and A 2 imposes the unit length constraint on z a . We can also insert the parameterization 



(2.3) into Hi and obtain the Lagrangain 



(A21) 



There is now a natural mean field theory of Cq + Ci which should yield all 4 phases of 
Ref. [26j The approximations are: 

• Ignore the gauge field A a r . 

• Factorize the 4-Fermi term, K± into N "^^^,^^ The field A^ a is to be determined 
self-consistently, and will be site-dependent. 

• Factorize Li into fermion and boson bilinears, as indicated by the parentheses. 

• Phases A and C will also have a <fi condensate. It should be sufficient to work with 
((f)) = in phases B and D, and determine their boundaries to phases A and C 

• Phase D should have N a = 0, and also (z* a z*p) = and (ify + i(;J) = 0. 

• The value of A 2 is determined as usual by solving the unit length constraint on z a . 



Appendix B: Square lattice antiferromagnetic in an applied gauge flux 



This appendix will carry out a computation similar to that of Section |II A| using gauge 



theoretical formulation in Eq. (2.5). Rather than a slowly varying Neel order n a (r) as 
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in Eq. (2.2), this appendix will have a slowly varying gauge potential A(r). The results 



here will be connected to those of Section II A via Eq. (2.7). However, a precise quantitative 



equivalence between Eqs. (2.2) and (2.5) requires inclusion of the last two terms of Eq. (A21 ) 



in Eq. (2.5), which we will not account for here. The importance of these omitted terms 



should be clear from Appendix A of Ref. |271 



Now we expand Eq. (2.5) to first order in Ay, and using Eq. (2.6) we can write H 



H + Hi where H has the same form as Eq. (2.9) but with the ip± fermions 



#o = ^(e k ^ t (k)^(k)+m^ t (k + Q) ( T^(k)) 



(Bl) 



while Hi in Eq. (2.11) is replaced by 



Hi = -i^t(r< - r^Aij (vl^j - ^}a z ^ 



E 

k,q 



A(q) 



dk 



^(k + q/2)<r*V(k - q/2) + 0(q 2 ) 



(B2) 



Note that Hi does not include the omitted terms represented by the ellipses in Eq. (2.5), 



which appear as the last two terms in Eq. (A21); this will be significant below. 



Now we will use the Kubo formula to determine the response to the applied gauge field 
in Hi. We will work to linear response order A, and to linear order in q. 

We have to carefully define an observable: it should be gauge invariant and spin-rotation 
invariant. For this reason we look at the response in the following 



My = ijje^^ijj 



(B3) 



We want to compute the change in (My) to linear order in A(q), and in the limit of small 
q. In momentum space 



(My) = ^e-^ + ^<^( k );^(p)) 
k, P 



(B4) 



r, - n) e 



ik r, r : <Q r I ,\ 



In the second term, we have assumed we are expanding to linear order in A, and so assumed 
momentum conservation in the fermion bilinear expectation value. 

We have to expand the first term to linear order in A, and so we expand the second term 
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in Eq. (B4) using Wick's theorem. 



(V/(k'W'(p) l V'(k + q/2)<r*V(k - q/2)> 



5 P ,k+Q+ q /25k',k- q / 2 G'(k - q/2)F(k + q/2) 

+ 5 P ,k+ q /25k',k+Q- q / 2 G'(k + q/2)F(k - q/2) 



Also from Eq. (2.14) 



<^(k)a^(k + Q)) = -2^F(k) 
Putting everything together 



6(Mij) = 2 e- ik ' (ri - r ^e iq ' (r ^ +ri)/2 e iQ -^A(q) • 

k, q ,w 



(B5) 



(B6) 



de k 
dk 



de 



G(k - q/2)F(k + q/2) + ^^(k + Q + q/2)F(k - q/2) + 



dF(k) 
dk 



(B7) 



Explicit evaluation shows that the expression in the square brackets does indeed vanish at 



q = 0, as is required by gauge invariance. Now expand Eq. (B7) to first order in q and find 

(B8) 



8{Mij) = 2^ e - ik - (ri - r ^e iq - (r ^ +ri)/ V Q ^A(q) -I(k,q) 

k, q 



where 
I(k,q) 



deu ( de 



dk 



q 



k+Q 



dk 



ds 



k+Q 



dk 



q 



dk 



E 



m/2 



-iu + E lk ) 2 (-iu + £ 2k ) 2 



(B9) 



Combining (B8) and (B9), we have the result analogous to Eq. (2.22): 



(k)c(k + Q)> = -2z.F(k) (d x A y - d y A a 



(BIO) 



where 



T-(k) 



de 



k+Q 



dk dk 



E 



m/4 



-z'u; + .Ei k ) 2 (-2'w + i^k) 



m / fek+Q <9e k \ (sgn(£ lk ) - sgn(£ 2k )) 



x 



dk dk J (£ lk - £ 2k ) 



(Bll) 



We have written Eq. (BIO) in terms of the original electron operators c(k): we are working 
to linear order in A, and so this order all variables can be mapped onto the original gauge- 
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invariant operators. Comparing Eq. (Bll) with Eq. (2.23), and using Eq. (2.7), we should 
expect equality between .F(k) and ^(k). However, while both functions have an identical 
symmetry structure, and similar singularities at possible Fermi surfaces (which is all we 
need), they are not precisely equal. This can be traced to the absence of precise equality 
between Eqs. (2.2) and (2.5), due to the omission of the last two terms in Eq. (A21), which 
were also important in previous computations^. 
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